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Although ab-initio calculations of relativistic Brueckner theory lead to large scalar isovector fields 
in nuclear matter, at present, successful versions of covariant density functional theory neglect the 
interactions in this channel. A new high precision density functional DD-ME5 is presented which 
includes four mesons a, to, 8, and p with density dependent meson-nucleon couplings. It is based to a 
large extent on microscopic ab-initio calculations in nuclear matter. Only four of its parameters are 
determined by adjusting to binding energies and charge radii of finite nuclei. The other parameters, 
in particular the density dependence of the meson-nucleon vertices, are adjusted to non-relativistic 
and relativistic Brueckner calculations of symmetric and asymmetric nuclear matter. The isovector 
effective mass m* — m* derived from relativistic Brueckner theory is used to determine the coupling 
strength of the 5-meson and its density dependence. 

PACS numbers: 21.10.Dr,21.30.Fe,21.60.De,21.60.Jz,21.65.Cd,21.65.-f,21.65.Ef 



I. INTRODUCTION 

Structure properties of nuclei described in the frame- 
work of effective mean-field interactions are remark ably 
successful over almost the entire periodic table fl2j] - 
Relativistic and non-relativistic versions of this approach 
enable an effective description of the nuclear many-body 
problem as an energy density functional. These energy 
functionals are usually adjusted to a variety of finite nu- 
clei and infinite nuclear matter properties. Although all 
these effective interactions are based on the mean-field 
approach, some differences will generally appear between 
them due to the specific ansatz of the density depen- 
dence adopted for each interaction. For instance, pre- 
dictions in the isovector channel of existing functionals 
differ widely from one another and, as a consequence, 
the density dependence of the symmetry energy is far 
from being fully determined. This has an impact on fi- 
nite nuclei properties as, for example, the neutron skin 
thickness. Mean field models which accurately describe 
the charge radius in 208 Pb, predict neutron radii between 
6.6 and 5.8 fin. For these reasons, one of the main goals in 
Nuclear Physics is to build a universal density functional 
theory based on microscopical calculations [lj, [l4| • This 
functional should be able to explain as many as possi- 
ble measured data within the same parameter set and 
to provide reliable predictions for properties of nuclei far 
from stability not yet or never accessible to experiments 
in the laboratory. It should be derived in a fully micro- 
scopic way from the interactions between bare nucleons. 
At present, however, attempts to derive such a density 
functional provide only qualitative results for two rea- 
sons: first the three-body term of the bare interaction is 



not known well enough and second the methods to derive 
such a functional are not precise enough to achieve the 
required accuracy. Note that a 1 %c error in the binding 
energy per particle of symmetric nuclear matter leads to 
an error of several McV in the binding energy of heavy 
nuclei, an error which is an order of magnitude larger 
than required by astrophysical applications. Therefore, 
at present the most successful functionals are derived ei- 
ther in a fully phenomenological way from a very large 
set of experimental data, as for instance more than 2000 
nuclear binding energies 0] or, more recently, by an ad- 
justment to a combination of microscopic results and to 
a set of characteristic experimental data fl5l - fl8l |. 

With the same goal, in the more recent past a different 
approach was tried. Baldo, Schuck and Vinas employed 
the same route as in condensed matter physics and con- 
structed a functional (BCP) [l2[ where the bulk part is 
based on the fully microscopic calculations of Ref. [19| . In 
this more fundamental calculations, Baldo and collabo- 
rators investigate the nuclear and infinite neutron matter 
on the basis of the Bethe-Brueckner approach including 
three-body correlations of the Bethe-Faddeev type (20j . 
Their results for the Equation of State (EoS) of sym- 
metric and neutron matter are believed to be among the 
most accurate in the literature. Then, the BCP func- 
tional took as a benchmark the calculations of Ref. 
by means of a polynomial fit. In this way an accurate 
and analytic EoS as a function of neutron and proton 
densities were constructed covering the whole range from 
symmetric nuclear matter to pure neutron matter in den- 
sity ranges from zero to about two times saturation den- 
sity. Subsequently a finite range surface term dependent 
on three parameters together with the strength of the 
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spin-orbit term (four parameters in total) were added to 
the bulk part of the functional. The pairing correlations 
needed to describe open-shell nuclei were accounted by a 
density dependent <5-force with an effective mass equal to 
the nucleon mass that simulates pairing calculations in 
symmetric nuclear matter computed with the Gogny in- 
~ in the T 
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1 channel. Adjusting these free 
parameters to some selected nuclear experimental data 
yielded excellent results for nuclear masses and reason- 
able charge radii of the whole nuclear chart. In addition 
it has also been shown that the deformation properties of 
BCP functionals are similar to the ones found using the 
Gogny D1S force [HI, HH in spite of the fact that both 
models are clearly different. A recent review on the BCP 
functionals can be found in Ref. [13] 

In general, symmetries of nature help to reduce the 
number of parameters and to simplify the description. 
One of the underlying symmetries of QCD is Lorentz in- 
variance and therefore covariant versions of density func- 
tionals are of particular interest in nuclear physics. This 
symmetry allows one to describe the spin-orbit coupling, 
which has an essential influence on the underlying shell 
structure in finite nuclei, in a consistent way. Moreover 
it also puts stringent restrictions on the number of pa- 
rameters in the corresponding functionals without reduc- 
ing the quality of the agreement with experimental data. 
Self-consistent mean-field calculations starting from rela- 
tivistic Lagrangians have been very successful in describ- 
ing nuclear properties [25l - [29j . They arise from a mi- 
croscopic treatment of the nuclear many-body problem 
in terms of nucleons and mesons carrying the effective 
interaction between nucleons. Moreover, since the the- 
ory is relativistically invariant and the field and nucleon 
equations of motion arc solved sclf-consistcntly, they pre- 
serve causality and provide a self-consistent description 
of the spin-orbit term of the nuclear effective force and 
of the bulk and surface parts of the interaction. In addi- 
tion these functionals include nuclear magnetism [30| . i.e. 
a consistent description of currents and time-odd fields, 
important for odd-mass nuclei |3lj . excitations with un- 
saturated spins, magnetic moments [H| and nuclear ro- 
tations (33d . No new parameters are required for the 
time-odd parts of the mean fields. In non-relativistic 
functionals the corresponding time-odd parts are usu- 
ally difficult to adjust to experimental data and even if 
there are additional constraints derived from Galilean in- 
variance and gauge symmetry [34| these constraints are 
usually not taken into account in the successful function- 
als commonly used in the literature. The earlier ver- 
sions of covariant densit y fu nctional theory were based 
on the Walccka model [3a - l38j with phenomenological 
non-linear meson-interactions proposed by Boguta and 
Bodmer [3i| introducingin this way a phenomenologi- 
cal density dependence [1, d, |4(j ■ Later the non-linear 
models have been replaced by an explicit density de- 
pendence of the meson-nucleon vertices. This density 
dependence has first been determined in a phenomeno- 
logical way @, 0, G3 These models have shown consid- 



erable improvements with respect to previous relativis- 
tic mean-field models in the description of asymmetric 
nuclear matter, neutron matter and nuclei far from the 
stability valley. On the other hand one has tried to de- 
rive this density dependence in a microscopic way from 
Brueckner calculations in nuclear matter at various den- 
sities [4lTj44l ]. An example is Density Dependent Rela- 
tivistic Hadron Field theory [iH where the specific den- 
sity dependence of the meson-nucleon vertices is mapped 
from Dirac-Brueckner calculations where the in-medium 
interaction is obtained from nucleon-nucleon potentials 
consistent with scattering experiments. Therefore, if this 
ansatz is adopted, the effective theory is derived fully 
from first-principle calculations. Of course, the accuracy 
of the results obtained in this way is by no means satis- 
factory for modern nuclear structure calculations and a 
fit of additional free parameters is still needed. This fact 
allows to constrain the different possibilities and keeps 
the compatibility, at least theoretically, with more fun- 
damental calculations of infinite nuclear matter. 

As mentioned, there exist ab initio calculations of the 
nuclear EoS over a wide range of nucleon densities, i.e. 
far from densities currently reachable at the laboratory. 
In this sense, apart from the experimental data needed 
in the fitting procedure for determining an effective in- 
teraction, further steps on building a universal density 
functional may need to implement such ab initio infor- 
mation as it was done in the BCP model. Therefore, this 
EoS calculated from first principles can be understood 
as a temporary benchmark at supra- and sub-saturation 
densities of the energy per particle at different asymme- 
tries. Furthermore, from a theoretical point of view, con- 
sistency is desirable between predictions of both theories. 
Regrettably, to make them compatible is not only a prob- 
lem of including the EoS derived from realistic nucleon- 
nucleon potentials within the fitting procedure. It is also 
a problem of taking into account the proper density de- 
pendence in the different terms of the functional. For 
that, first-principle calculations are also thought to be 
the best candidate to help in building a universal en- 
ergy density functional, at least the bulk part of it since 
many-nucleon calculations are not feasible yet. Hence, 
for an effective and self-consistent treatment of the nu- 
clear many body problem, we propose here a novel and 
improved relativistic mean-field interaction with an ex- 
plicit density dependence of the meson-nucleon vertices 
in all the four spin-isospin channels compatible with fully 
microscopical calculations. 

The essential breakthrough of density functional meth- 
ods in the description of quantum mechanical many-body 
problems was Kohn-Sham theory [45|, [46| , where the ex- 
act density functional E[p] of Hohenberg and Kohn [47j 
was mapped in an exact way on an effective potential V c s 
in a single particle Schroedinger equation, which forms 
the starting point of all modern applications of density 
functional theory. In covariant density functional theory 
this effective potential corresponds to the self-energy in 
the Dirac equation, which can be decomposed into four 
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channels characterized by the relativistic quantum num- 
bers of spin and isospin, the scalar isoscalar channel (J = 
0, T = 0), the vector isoscalar channel (J = 1, T = 0) the 
scalar isovector channel (J = 0, T = 1) and the vector 
isovector channel (J = 1,T = 1). In the Walecka model 
these channels arc connected with the exchange of mesons 
carrying the corresponding quantum numbers. 

However, in nearly all the present successful phe- 
nomenological applications of covariant density func- 
tional theory to nuclear structure based on the relativis- 
tic Hartree model only the three mesons a , u) , and 
p are taken into account. The scalar isovector <5-meson 
(ao (980)) causing a splitting of the effective mass between 
protons and neutrons is neglected, because it has turned 
out that usual data such as binding energies and radii of 
finite nuclei do not allow to distinguish scalar and vector 
fields in the isovector channels. Allowing independent pa- 
rameters for the p- and the J-mesons leads to redundan- 
cies in the fit. By the same reasons also modern relativis- 
tic Hartree-Fock [48l - l50T ] and Hartree-Fock-Bogoliubov 
[5~J ] calculations neglect the (5-meson in the Lagrangian 
and, as a consequence, in the direct term. Of course 
the Fock term of these calculations contains also con- 
tributions to the scalar isovector channel. Microscopic 
investigations by Huber efc al [52 . |53j and phenomeno- 
logical studies |5J-[64| in the the literature stressed that 
mean-field models which neglect the (5-meson arc likely 
to miss important ingredients in describing properly very 
asymmetric nuclear matter, in particular at high densi- 
ties. The proton fraction of /3-stable matter in neutron 
stars can increase and the splitting of the effective mass 
can affect transport properties in neutron stars and heavy 
ion reactions. However, as long as the parameters of this 
meson are not fixed, such investigations are somewhat 
academic. Therefore we derive in this manuscript the 
<5-nucleon vertex and its density dependence from mod- 
ern microscopic calculations based on the bare nucleon- 
nucleon force of the Tubingen group [65j . 

The relativistic mean field model DD-ME<5 obtained in 
this way is an extension of the DD-ME model developed 
by the Munich group 0, [H3] based on the density depen- 
dent relativistic Hartree theory. The DD-ME model has 
the following degrees of freedom: the proton, the neu- 
tron and three mesons carrying the nuclear interaction, 
namely a-, uj- and p-mesons. In addition to these degrees 
of freedom, we include here a new one, the (5-meson by 
the reasons pointed out before. Since this meson provides 
a treatment of the isospin more close to the microscopic 
investigations we can hope that it improves the reliability 
of the models for predictions in nuclei far from stability 
with large isospin - planned to be studied experimentally 
at the new Rare Ion Beam Facilities [66[. Apart from 
the inclusion of the (5-meson in DD-ME<5 the DD-ME<5 
model differs from the earlier DD-ME models in that the 
parameters of DD-ME were all adjusted to experimental 
data based on finite nuclei properties, whereas those of 
DD-ME(5 are largely based on microscopic ab-initio cal- 
culations in nuclear matter. Only four of the parameters 



of DD-ME<5 are fitted to finite nuclei. 

The paper is organized as follows: after establishing 
the DD-MEJ model in Sect. HH we discuss in Sect. UTE] 
the strategies to determine the parameters of the La- 
grangian and compare in Sect. HVl the results of this novel 
effective interaction with the experiment and with the 
non-rclativistic BCP model [l2| and the completely phe- 
nomenological DD-ME2 model , in particular binding 
energies, charge radii and neutron skins in spherical nu- 
clei. 



II. DENSITY DEPENDENT HADRON FIELD 
THEORY 

A. Lagrangian and equations of motion 

Density dependent relativistic hadron field theory 
which forms the basis of the DD-MEi5 interaction has 
been formulated and extensively discussed in Refs. [42|, 
l67l . Hll . Here we present only the essential features of the 
mean-field equations of motion derived from such a the- 
ory. The relativistic Lagrangian includes neutrons and 
protons represented by the Dirac spinors ip of the nu- 
cleon, the four mesons (cr, u, S and p) carrying the effec- 
tive nuclear strong interaction represented by the fields 
a, <5, and fp*, and the photon field accounting for 
the electromagnetic interaction. The index p indicates 
the time- and space-like components of the vector fields 
and the arrow indicates the vector nature of a field in 
isospin space. As mentioned, the (5-meson should be in- 
cluded if one wants to follow the theoretical indications of 
Dirac-Brueckner calculations in asymmetric nuclear mat- 
ter and so we do. The Lagrangian has the following parts 

C = C N + Cm + £int (1) 

where C N is the nucleonic free Lagrangian 

C N = i> (*7^ M - m) ^ (2) 
Cm is the Lagrangian of free mesons 

Cm = \ (d^aWa - m> 2 ) + ^ (dJd^S- m» (3) 

- ^ V 

and Cm is the Lagrangian describing the interactions. 
Its algebraic expression is 

C int = gaTpvip + gs^rSijj 

where m is the nucleon mass (commonly taken as 939 
McV), the field strength tensors for the vector fields are 

n"" = d»u v - d"^, (4) 
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and correspondingly R^ v and . The electric charge is 
e for protons and zero for neutrons. The meson-nucleon 
vertices are denoted by <jj for i = a, w, S, and p. Since 
covariance is required and the quantity y^jpj^ is in the 
rest frame identical to the baryon density p = p n +p p , the 
nucleon-meson vertices generally depend on this quantity. 
Because of the relatively small velocities the difference 
between ^[Jp~j^ an d P 1S negligible in all practical appli- 
cations. The subindcx n or p is used to indicate whether 
we are considering neutrons or protons respectively. 

The equations of motion are derived from the classi- 
cal variational principle and we obtain for the nucleon 
spinors the Dirac equation 



state expectation values defined as, 



['T M (i& t -J^)-m*]^ = 



(5) 



where m* = m — S s is the effective Dirac nucleon mass 
and S M and S s are the vector and scalar self-energies 
defined as follows, 



£ s = 9<j{p) g + 9s(p)tS 



(6) 
{7) 



here, (0) indicates the usual definition of the vector self- 
energy and (r) the rearrangement term of the vector self- 
energy 



E<°>" = g u (p)w» 
v (r) _ dg, 



9 P (p)T 3 P%+eA» 



(8) 



dp 



ap 



dgu o d 9s , s . dg p 



Here eA° is the direct term of the Coulomb potential. 
As in most RMF models we neglect in these investiga- 
tions the Coulomb exchange term which plays an impor- 
tant role in pn-KPA calculations [7(|. The static mean 
field approximation used throughout this investigation 
preserves the third component of the isospin. As a con- 
sequence the other two components of the densities and 
fields carrying isospin vanish. In Eq. [5] and the following 
equations p!] represents the time-like component of the 
p- meson field, whereas p 3 = p n — p p and p% = p a n — p p rep- 
resent the isovector part of the baryon density and of the 
scalar density. The rearrangement term is a contribution 
to the vector self-energy due to the density dependence 
of the meson-nucleon vertices. The equations of motion 
for the mesons arc, 



{d v d v 

{p v d v + 
{p v d u - 
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P% = 



~g*(p) p s 

+9u{p) f 
-9s(p)pl 
+9 P {p) 3a 
+e Jit 



(10) 

(11) 
(12) 

(13) 
(14) 



d v d v A^ = 

where the different densities and currents are the ground 



p s = (o\W\o) 
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= Pn 


+ p s P 


(15) 


f = (oi^yvio) 


= 7 M 
Jn 
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pI = (oiv^io) 
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= Pn 


- Pp 


(17) 


% = (0|^r 3 V|0) 


Jn 




(18) 



B. Asymmetric infinite nuclear matter 

1. Energy density and pressure 

In infinite nuclear matter we neglect the electromag- 
netic field. Because of translational invariance, the Dirac 
equations can be solved analytically in momentum s pac e 
and we obtain the usual plane- wave Dirac spinors [69(. 
Filling up to the Fermi momenta kp T for r = n or p, we 
find the densities 
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(2tt) 3 J\ M< k PT E T (k) 



m T 
2^2 



kp T Ep T — m* 2 In 
and the meson fields 



(19) 
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+ Pp) 


(21) 
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ml 


-P S P ) 


(23) 


9 P (p) { 
ml [Pn 


-P P ) 


(24) 



where E T (k) = -\/k 2 + m* 2 and where the Fermi energy 
of neutrons and protons is given by Ep T = E T (kp T ). 
Now, we calculate the energy density (e) and pressure 
(P) from the energy-momentum tensor, 



rpp.1. 



dC 



d u fa - trc 



d{dn<f>i) 

where fa runs over all possible fields, 
e = (0|r 00 |0) 

J [3E Fn p n + m* n p s n ] + j [3E FpPp + m*p s p ] 



(25) 
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and 
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1 3 



|T"|0) 



(27) 



[ E FnPn ~ Tn*nPn 



[E Fp p p - m* v p s p 



[mW-mK^f+mlSl-mKpVf] 



+ (p„ + P P )£ (r)0 

Only the pressure has a rearrangement contribution. We 
have checked that the pressure derived from the energy- 
momentum tensor coincides with the thermodynami- 
cal definition: p = p 2 [d(e/ p) / dp] and that the energy- 
momentum tensor is conserved dpT^ = 0. 



The quantity u a (p,m*) will be needed below. In these 
equations we used the fact that for a = we have p n = 
p p = p/2, p s n = p s p = p s /2, p% = 0, m* = m* p and 

we have defined kp = 3ir 2 p/2 and Ep = J k 2 F + m* 2 . 

In symmetric nuclear matter, the effective mass and the 
scalar density read 



* 9t7 S 

m = m Tr-p 

mi 



*2 , i k F + E F 
kfEf — m In 



(35) 
(36) 



respectively. Close to the saturation density us ~ 1 is 
a very good approximation and we find in this case an 
analytical approximation for (p) 



2. The symmetry energy: S2(p) 

Assuming charge symmetry for the strong interaction 
(the nn and pp interactions are identical but different, 
in general, from the np interaction) , the total energy per 
particle in asymmetric nuclear matter can be written as 
follows, 

- = ^ =e(p,a) =e(p,a = 0) + S 2 {p)a 2 + (D[a 4 } (28) 
p A 

where p = p„ + p p is the baryon density and a = 
(Pn — Pp)/{Pn + Pp) measures the neutron excess. The 
term proportional to a 2 is the so called symmetry en- 
ergy of infinite matter and terms proportional to a 4 (and 
higher) can be neglected to very good approximation. 
The symmetry energy S2G0) is defined as, 



nS / \ 1 9 S ( m * 



(37) 



and, therefore, the contribution to the symmetry energy 
coming from the nuclear strong interaction (potential 
part) as described by this kind of models can be writ- 
ten in the simple form, 



sr\p)=s^p)+s s 2 ( P ) 
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(38) 



The symmetry energy is often expanded around the sat- 
uration density p sat 



S 2 (p) 



1 f d 2 e(p,a) 

2 { da 2 



(29) S 2 (p) = J+- (p-p sat ) + 



a=0 



Models including the (5-meson provide a richer descrip- 
tion of the isovector sector of the nuclear strong interac- 
tion. For that reason it is important to understand its 
effects on asymmetric nuclear matter and for that we give 
the analytic expressions for the symmetry energy of the 
model discussed in the last section 56, 571: 



S 2 (p)=S? n (p) + S£(p) + Sfo) 



with 
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where for i = a, 8 
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where J is the symmetry energy at saturation, and L and 
K sym arc proportional, respectively, to the slope and the 
curvature of the symmetry energy at saturation. 
Using the analytical expressions (|3"3"|) we find 
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dp 
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where the functions Uj, v and w depend on p and m* 



v(p,m*) 
w(p,m*) 



y m* 
p dm* 
m* dp 

mi 



E F 



,£L_P_dga_ _P_ 

m* g a dp Ef 



(44) 
(45) 



The strength of the a- and w-nucleon vertices is quite 
well determined by experimental data as compared with 
the strength of the isovector meson-nucleon vertices. On 
the other side, with only the p-nucleon vertex, one is 
able to reproduce properties of finite nuclei [7l[ and to 
account for the symmetry energy around saturation in 
rather good agreement with available empirical indica- 
tions. However, to reproduce nucleon-nucleon scattering 
measurements in the vacuum, one needs to incorporate 
a scalar-isovector meson into the parameterization of the 
two-body nuclear interaction [72j ]. Microscopic deriva- 
tions of the nuclear fields using relativistic Brucckner 
theory [42l l52l l65l I73l - l75"j or non- relativistic Brueck- 
ner theory [43|, [44J show clearly that the scalar field in 
the nuclear interior has an isovector part. These reasons 
motivate one to incorporate the i5-meson also in models 
of covariant density functional theory and to study its 
influence on properties such as the symmetry energy, the 
effective mass splitting between protons and neutrons in 
asymmetric matter, the isospin dependence of the spin- 
orbit potential and the spin-orbit splittings far from sta- 
bility. 



Density dependence of the meson-nucleon 
vertices 



where p sa t is the saturation density of symmetric nuclear 
matter and x = p/p s &t- For the functions fi(x) we follow 
Refs. [f| 0, E3] and use the Typel-Wolter ansatz: 



l + bi{x + dif 
1 + Ci(x + ei) 2 



(47) 



As in Refs. @, E3 El we use the value p sat = 0.152 
fm -3 . In fact, this choice is very close to the saturation 
density obtained in the following fit. As we see from the 
ansatz (|47p the actual value of p sat is irrelevant for the 
calculations. It can be completely absorbed in the values 
of the parameters bi , c,; , di , and ej . It is only used to make 
them dimensionless. By definition, the parameters at are 
constrained by the condition /(l) = 1. In the earlier 
applications [g, 0, Ell this ansatz was only used for the 
a- and w-meson. The density dependence of the isovector 
coupling g p (p) was described by an exponential and the 
(5-meson was neglected. Here we use the same ansatz (|4"7|) 
also for the isovector mesons S and p. This turned out to 
be necessary in order to obtain a density dependence of 
the 5-coupling similar to that derived from microscopic 
ab-initio calculations in Refs. [42|, EH E3 • We impose as 
in Ref. @ the constraints e a = d a , e w = d u , f a [x = 1) = 
f u (x = 1) and f { (x = 0) = 0. We work with meson 
masses m u = 783 MeV, mj = 983 MeV and m p = 763 
MeV. The nucleon mass is m = 939 MeV. All in all, 
the model has 14 adjustable parameters. Namely, the 4 
coupling constants gi(p sa ,t) in the 4 relativistic channels 
(Lorentz-scalar, Lorentz- vector, isoscalar and isovector), 
9 parameters describing the density dependence in the 
functions fi(x), and the cr-mass m a allowing for a finite 
range and a proper description of the nuclear surface. 



Here we describe the density dependence of the meson- 
nucleon vertices used for the new interaction DD-ME<5. 
We start from modern fully microscopic calculations in 
symmetric nuclear matter and pure neutron matter at 
various densities and try to determine the density depen- 
dence of the vertices gt(p) by fitting to those data. Of 
course, it is well known that successful density function- 
als can, at present, not be determined completely from 
ab-initio calculations. Therefore, we introduce in the fit 
not only results of microscopic calculations but also a set 
of data on binding energies and radii in specific finite 
nuclei. 

In a first step we have to choose a form of the den- 
sity dependence of the various vertices, which is flexi- 
ble enough to reproduce the microscopic calculations. In 
Refs. [U EH the meson-nucleon vertices of density de- 
pendent RMF theory have been related to the scalar and 
vector self-energies obtained from Dirac-Brucckncr (DB) 
calculations in infinite nuclear matter. The density de- 
pendence deduced from DB calculations is 

9i(p) = 9i{Psa,t)fi(x) for i = a,Lj,8,p (46) 



D. Calculation of finite nuclei 

The self-consistent results for masses include a micro- 
scopic estimate for the center-of-mass correction: 

(P 2 ) 

E cm = - { -^ (48) 
2mA 

where P cm is the total momentum of a nucleus with A 
nucleons. The expression 

r c = ^(r 2 }+0.& 2 (49) 

is used for the charge radius. The description of open 
shell nuclei requires pairing correlations. We introduce 
this through the BCS approach with a seniority zero force 
in the soft pairing window described in Ref. [78j |. For 
the fit of the parameters of the Lagrangian described in 
the next section the constant gap approximation (79j has 
been used and the gap parameters have been derived from 
the experimental binding energies by a 3-point formula. 
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III. THE PARAMETERS OF THE 
FUNCTIONAL DD-ME5 

In this section we describe the determination of the 
parameters of DD-ME<5. Earlier fits of rclativistic La- 
grangians have shown that the usual set of experimental 
ground state properties in finite nuclei, such as binding 
energies and radii do not allow to determine more than 
7 or 8 parameters 0. Two of them (g a / m a and gu/niu) 
determine the saturation energy and the saturation den- 
sity of symmetric nuclear matter [37| . one of them (m a ) 
is fixed by the radii in finite nuclei and another one of 
them (g p /m p ) determines the symmetry energy J at sat- 
uration. The additional parameters (as for instance g2 
and 93 in the- nonlinear meson coupling models NL1 3] or 
NL3 [J] or the three parameters in the ansatz (|47| for den- 
sity dependence in the isoscalar channel of DD-ME1 Q 
and DD-ME2 are determined by the isoscalar sur- 
face properties and are necessary to describe deforma- 
tions and the nuclear incompressibility properly. Finally 
one parameter (a p in DD-ME1 or DD-ME2) is needed to 
describe the density dependence of the symmetry energy 
by a fit to the experimental data on the neutron skin 
thickness. 
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FIG. 1: (Color online) Equation of state for symmetric nu- 
clear matter and for pure neutron matter as a function of the 
nucleon density. The dots represent the predictions of the 
BHF calculations and the line our fit to reproduce this 
data. 

In order to calibrate the 14 free parameters of the DD- 
ME5 functional we therefore added pseudo-data in the 
form of results of modern microscopic non-rclativistic and 
rclativistic Brueckncr calculations. To this end, we se- 
lected the EoS of symmetric nuclear matter and of neu- 
tron matter (sec Fig. [1} derived by Baldo ct al. [l!| in 
a state-of-the-art non-relativistic Brueckner calculation 
including relativistic corrections and three-body forces. 
We also used as a benchmark the isovector part of the 
effective Dirac mass m* — to* (see Fig. [2| derived by the 
Tubingen group (6o| in relativistic Dirac-Brueckner the- 



ory. The use of non-relativistic results for the EoS and 
of relativistic results for the isovector effective mass may 
seem somewhat arbitrary. However we have to keep in 
mind that the non-relativistic calculations of the Cata- 
nia group are more sophisticated than presently avail- 
able Dirac-Brueckner calculations, because they include 
not only relativistic effects but also three-body forces. 
On the other side it is very complicated to deduce Dirac 
masses from a non-relativistic calculation which does not 
distinguish between Lorentz scalars and vectors. This is 
in principle possible [43|, [44| , but it is difficult and con- 
nected with additional uncertainties. With this caveat 
in mind, we decided to use a reliable relativistic Brueck- 
ner calculation [6^] providing directly the effective Dirac 
masses to* and to* and their difference, a quantity di- 
rectly connected with the scalar isovector part of the self 
energy. The isovector part of the effective Dirac mass 
to* — m* depends only on the 5-meson. It vanishes for 
all the conventional Lagrangians without <5-meson. The 
density dependence of this quantity is therefore the opti- 
mal tool to get information about the density dependence 
of the S- meson vertex gs(p). 




FIG. 2: (Color online) Proton-neutron effective mass splitting 
as a function of the nucleon density in pure neutron mat- 
ter. The dots represent the predictions of the DBHF calcula- 
tions [HHl and the line our fit to reproduce this data. 



Keeping this in mind, we determine 10 of the 14 param- 
eters in the Lagrangian of DD-ME<5 by these pseudo-data 
obtained from ab-initio calculations of nuclear matter. 
These parameters define the density dependence for the 
various meson-nucleon vertices (9 parameters) and the 
strength gs(p s &t) of the 5-meson. Only a reduced set of 
4 parameters (.^(psat), ff w (/O sa t), g P {p sa ,t), and m a ) are 
fitted to the masses and charge radii of finite nuclei. 

A. Strategy of the parameter fit 

Since the mean field equations of motion have to be 
solved self-consistently, we need a good starting param- 
eter set before fixing the 5-mcson coupling to the above 
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mentioned calculations fully. The density-dependent me- 
son coupling model DD-ME2 provides us with an 
excellent description of nuclei all over the periodic table. 
Though DD-ME2 neglects the 5-meson, it is based in the 
isoscalar channel on the same ansatz ([47]). Therefore, we 
used DD-ME2 as a starting point of our investigations. 
We proceeded in three steps: 

1) In the first step we performed an overall fit with all 
the 14 parameters. For the data we have chosen on 
one side the three microscopic curves for the EoS 
in Figs. [1] and [2] and on the other side the same 
set of data of finite nuclei which has been used in 
Ref. [Io| for the determination of the parameter 
set DD-ME2 (see Table II of this reference), i.e. 12 
binding energies of spherical nuclei distributed all 
over the periodic table and 9 charge radii. Due to 
the fact that the density dependence in the isovec- 
tor channel is determined by the equation of state 
of neutrons it was not necessary to include neutron 
skin thicknesses (r„ — r p ) data. This fit provides 
us with a relatively stable starting point for a sub- 
sequent fine tuning of the model. Moreover, since 
the J-meson is little influenced by the overall fit to 
finite nuclei both gs(p sa ,t) and fs(x) (4 parameters) 
are relatively well determined already in this step 
and we need only a fine tuning of the remaining 
parameters in the next two steps. 

2) In the second step we keep the four meson masses 
and the four parameters describing the density de- 
pendent vertex gs(p) of the 5- meson fixed and de- 
termine the 9 parameters describing the Typel- 
Woltcr ansatz for the density dependent vertices 
gi(p) of the remaining three mesons (i = a, u> and 
p) by a very accurate fit to the nuclear matter data 
shown in Figs. [T] and [2j Involving only nuclear mat- 
ter data, this is a relatively fast calculation and 
as a result we obtain the three density dependent 
vertices g%{p) for i = a, u, p. In this way we de- 
scribe with high precision the EoS for symmetric 
nuclear matter and pure neutron matter as well 
as the isovector part of the effective Dirac mass 
Am* = m* — m* . 

3) In the last step we keep the ^-meson parameters 
as determined in step 1 and the density-dependent 
functions fi(x) for i = a, w, p are frozen at the val- 
ues found in step 2. We refine the remaining 4 
parameters g a {p sa t), 9u(fiaat), 9p(p a et), and m a to 
the binding energies of 161 spherical nuclei and the 
charge radii of 86 nuclei shown in Table [V] taking 
into account in this fit also the pseudo-data of the 
nuclear matter properties used in step 1 and 2 with 
certain weights. It turns out that the values of 
.9<r(Psat), ff w (Psat), g P (p sa ,t), and m CT obtained in this 
fit differ only slightly from the values determined in 
step 1 and that the nuclear matter results (EoS in 
symmetric nuclear matter and pure neutron mat- 
ter as well as the isovector Dirac mass) differ only 



marginally from the results obtained in step 2. As 
a consequence the procedure involving step 2 and 
3 does not have to be repeated. 



i rrii (MeV) ffi(p sa t) a% b t a d t e. 



a 566.1577 10.3325 1.3927 0.1901 0.3679 0.9519 0.9519 

w 783.0000 12.2904 1.4089 0.1698 0.3429 0.9860 0.9860 

5 983.0000 7.1520 1.5178 0.3262 0.6041 0.4257 0.5885 

p 763.0000 6.3128 1.8877 0.0651 0.3469 0.9417 0.9737 



TABLE I: The parameter set DD-ME5 with the 5-meson. 
It includes 14 independent parameters. Only 4 of them 
(g<r(psat), g^(psat), Sp(paat), and ma) are fitted to finite nu- 
clei. The other 10 are derived from an adjustment to ab-initio 
calculations in infinite nuclear matter [ltj, [65(] 

The final parameter set DD-ME<5 obtained in this way 
is given in Table HI It is compared with the parameter 
set DD-ME2 in Table HU We observe a large difference 
in the value of the the p-nucleon vertex g p . This can be 
understood by the fact that we have in DD-ME2 only one 
meson (p) in the isovector channel with a very different 
density dependence. We also observe a very different 
density dependence for the isoscalar mesons a and u. 
All in all this has only a minor effect in the low density 
region p < p sat but it has a large effect at high densities 
as shown in Fig. [3] 



i rrii (MeV) #i(p sa t) <H b, a di e. 



a 550.1238 10.5396 1.3881 1.0943 1.7057 0.4421 0.4421 
lj 783.0000 13.0189 1.3892 0.9240 1.4620 0.4775 0.4775 
p 763.0000 3.6836 0.5647 



TABLE II: The parameter set DD-ME2 which does not con- 
tain a 5-meson. It includes 8 independent parameters. All are 
fitted to finite nuclei. In DD-ME2 the density dependence of 
the p-meson is given by g p {p) = 3p(p sa t) exp(— a p (x — 1)) 



B. \ 2 definition 

Our fit is performed through a \ 2 test of the form 

, "data 

x * = -i- £ w \ (0?° del ~ Off (50) 

2 — 1 

where nd ata is the number of data points and Wi the 
weight associated to each data point. G\ ei is the exper- 
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imcntal value for finite nuclei and the pseudo data ob- 
tained by ab-initio calculations in nuclear matter. The 
observablcs in finite nuclei used for the fit are the bind- 
ing energies of 161 nuclei and the charge radii of 86 nu- 
clei given in Table [V] All of the isotopes are spherical 
even-even nuclei and the data are taken from the liter- 
ature [goj]. In the standard definition of a x 2 test, the 
weights Wi should be inversely proportional to the exper- 
imental uncertainties. However, in the case of energies 
these are usually so small, that they cannot be used as 
relevant quantities. We therefore used the weights given 
in Table HTH i.e. Wi = 1/0.5 MeV -1 for the masses and 
Wi = 1/0.01 fm -1 for the radii. For the fit to the results 
of ab-initio calculations in nuclear matter we use ?idata 
mesh points in a certain density range (see Table ITTT|) and 
we assume a relative accuracy of 3 %. The minimization 
of x 2 is carried out by means of a variable metric method 
algorithm included in the MINUIT package of Ref. 



81 



a 


Wi 




p (fm 


- 3 ) 


x 2 


Ref. 


B 


1/0.50 MeV" 1 


161 






23.40 


[80] 


r c 


1/0.01 fin" 1 


86 






2.90 


[82] 


e(p,a = 1) 


1/(0.03 x e>o 


30 


0.01 - 


0.30 


3.42 


[19] 


e(p,a = 0) 


1/(0.03 x Oi) 


30 


0.01 - 


0.30 


7.03 


[19] 


m* - m* n 


1/(0.03 x Oi) 


25 


0.04- 


0.20 


0.39 


[65J 



TABLE III: Specifications of the \ 2 definition. The total x 2 
found for DD-ME5 is almost 40, the partial contributions to 
it are listed in the fifth column. 



In the first step of the fit discussed in the last section 
we minimize the quantity 



X(i) 



xl 



2,2 
Xrv ' Xsv 



2 

Xneut 



(51) 



At this stage all the 14 parameters of the model are varied 
and the data of finite nuclei are restricted to the masses 
and charge radii of the 12 nuclei used in the fit of the pa- 
rameter set DD-ME2 in Ref. [lo|. As wc have described 
in the last section, the parameters for the <5-meson ob- 
tained from this fit are no longer changed. We given 
them in the third line of Tabic fl] 
In the second step we minimize 



X(2) 



x: 



Xncut 



XAm* 



(52) 



for nuclear matter data with respect to the 6 constants 
characterizing the density dependence of the a, uj and p 
mesons and the 3 couplings g a {p sa t), 5w(Psat), g P {p sa t)- 
(Here, the J-meson and m a are hold at the values found in 
the first step.) The obtained values for the ai,bi,Ci,di, 
constants that define the density dependence of the i = 
u,uj,p meson- nucleon vertices are given in Table [IJ 



In the third step we minimize \ 2 f° r the nuclear matter 
data and the 161 binding energies and 86 charge radii 
given in Table |Vj 



*(3) 



xl 



\ r 



X 



sym 



v 2 

A ncut 



(53) 



Now we fit only a restricted set of 4 parameters, i.e., 
the 3 couplings g a (p sa t), g^(/O sa t), 5 P (p sa t), and m a . The 
resulting values are given in the first two columns of Ta- 
ble m 

We have to emphasize that only four free parameters 
.9<r(Psat), gu(p sa ,t), g P (p sa t) and m a have been used in the 
final fit to the experimental data in finite nuclei. The 
other 10 parameters are derived from ab-initio calcula- 
tions. This is in contrast to the typical rclativistic and 
non-relativistic fits of mean-field interactions, where com- 
monly around 10 free parameters are adjusted to data in 
finite nuclei. It is also worth to remember that adding 
the (5-meson has improved our theoretical picture of the 
nucleus and of the EoS of asymmetric nuclear matter. 

So far we have used in the fit only nuclei with spher- 
ical shapes. The pairing correlations are treated in the 
first step in the constant gap approximation with gap 
parameters derived from the odd-even mass differences. 
For a full description of nuclei all over the periodic table 
which includes also regions where the experimental bind- 
ing energies are not known, we introduce a more general 
description of pairing by means of a monopole force with 
a constant matrix element fitted to reproduce the exper- 
imental binding energies of the nuclei in Table [Vj Wc 
obtain for the set DD-MEJ the values G n = 32.44/A 
MeV for neutrons and G p = 29.76/A MeV for protons. 
In order to have a fair comparison for the results in finite 
nuclei we treated in the following the pairing properties 
of the set DD-ME2 also by a monopole force. In a sim- 
ilar fit we found for DD-ME2 the strength parameters 
G n = 29.86/^1 MeV and G p = 28.92/^1 MeV. In all these 
calculations the soft pairing window described in Ref. [78| 
has been used. 



IV. RESULTS 

A. Nuclear and Neutron Matter Equations of State 

The nuclear matter properties at saturation computed 
with the DD-ME(5 functional are given in Table IIVI 
These properties do not fully coincide with the ones of 
the fully microscopic calculation in [l|| . The reason for 
that is that in the microscopic calculation, the EoS is 
very flat around saturation density and some deviation 
between the microscopic results and the DD-ME<5 fit ap- 
pear. These differences remain within the uncertainty of 
the state of the art of present numerical microscopic cal- 
culations. They are too small to be seen on the scale of 
Fig. [TJ They are, however, important for a fine tuning of 
the results. 

In order to investigate the quality of the predictions of 
the density functional DD-ME<5 in the high density do- 
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main, we show in Fig.[3]the pressure (|27[) computed with 
this functional function of the density. It is com- 
pared with the pressure derived from the microscopic 
calculation of Ref. [l9| as well as with the results de- 
rived from the non-linear meson coupling; models NL3 [|| 
and FSUGold @ and from DD-ME2 M- We see that 
both, microscopic and DD-ME5 calculations, are within 
the shaded area which corresponds to the "experimen- 
tal region" estimated from simulations of heavy-ion col- 
lisions dH]. The standard non-linear a-uj model NL3 is 
outside of this region while the FSUGold model - which 
has an additional non-linear u-p coupling that softens the 
symmetry energy (see below) - is inside the shaded area 
in rather good agreement with DD-ME<5 and the micro- 
scopic results. The results of the parameter set DD-ME2 
are slightly outside of the shaded area. 

An important quantity in nuclear physics and astro- 
physics, directly related with the EoS of asymmetric nu- 
clear matter, is the symmetry energy (|30|) . The value of 
the symmetry energy derived from successful mean field 
models lies roughly in a window of 30 — 35 MeV at satura- 
tion. However, the density dependence of the symmetry 
energy is much more uncertain. This fact entails im- 
portant consequences for a number of isospin-dependent 
observables. As a paradigmatic example, one may recall 
that different accurate mean-field models which repro- 
duce well the binding energy and charge radius of the 
nucleus 208 Pb predict largely different values for the neu- 
tron skin thickness of this isotope, ranging from 0.1 — 0.3 
fm. This fact points out that the isovector properties of 
the different models are, actually, not well constrained 
by the binding energies and charge radii of stable finite 
nuclei used to fit the effective interactions. 

In nuclear mean field models, a strong linear correla- 
tion exists [m IIH between the size of the neutron skin 
thickness of a heavy neutron-rich nucleus such as 208 Pb 
and the L parameter defined in Eq. (|39[) . i.e., the slope of 
the symmetry energy at saturation. Recent constraints 
on the L parameter have been obtained using a variety 
of observables such as, for instance, isospin diffusion [86- 
HH and isoscaling [82-[93j in heavy ion reactions, some 
collective excitations in nuclei I7fi, 0^99 and the neu- 





DD-ME5 


DD-ME2 
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0.152 


0.152 


[fm- 3 ] 


€sat 


-16.12 


-16.14 


[MeV] 


Koo 


219.1 


250.89 


[MeV] 


J 


32.35 


32.30 


[MeV] 


L 


52.85 


51.26 


[MeV] 


m* jm 


0.609 


0.572 





TABLE IV: Nuclear saturation properties as predicted by the 
parameter sets DD-ME<5 and DD-ME2 
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FIG. 3: (Color online) The pressure defined in Eq. (f27|) for 
symmetric nuclear matter (upper panel) and for neutron 
matter (lower panel) as a function of the density. The re- 
sults for DD-ME5 are compared with those of the density 
functional NL3 from Ref. Q, FSUGold from Ref. 0, DD- 
ME2 from Ref. [lj}, and microscopic BHF calculations from 
Ref. [Tg| [ . The shaded area represents experimental results 
from Ref. El. 



tron skin thickness in finite nuclei [97J, |98j measured in 
antiprotonic atoms [9^, llOOj ] . The analysis of all these re- 
sults suggests that the L parameter is roughly within the 
window 45 - 75 MeV [97| ■ The new experimental efforts 
to measure the neutron radius of 208 Pb may turn out to 
be helpful to deduce in the future narrower constraints 
on the slope L of the symmetry energy through th e cor - 
relation of L with the neutron skin thickness ^a, 1 10 ill . 



The L value predicted by our novel DD-ME6 functional 
is 53 MeV, close to the resul t of the microscopic calcula- 
tion of 66.5 MeV in Ref. f 103| ] . The density dependence of 
the symmetry energy exhibited by DD-ME<5, is displayed 
in Fig. |H We see that DD-ME<5 predicts a rather soft 
density dependence of the symmetry energy which lies 
inside the shaded region derived from the empirical law 
Slip) = 31.6(p//9o) 7 MeV imposing the range discussed 
above: 45 MeV < L < 75 MeV. It turns out that the 
density dependence of DD-ME(5 and DD-ME2 is prac- 
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ME<5 and an exponential density dependence for DD- 
ME2 (see Eq. (7) in Ref. 0). 



0.2 0.3 



FIG. 4: (Color online) The symmetry energy in Eq. (|29|l as 
a function of the density. The results for DD-ME5 are com- 
pared with th ose o f the density functionals NL3 from Ref. , 
G2 from Ref. [To^], FSUGold from Ref. H, and microscopic 
DBHF calculations from Ref. [104| . The shaded area rep- 
resents the empirical region suggested by the available con- 
straints on the L parameter discussed in [97| . 



tically the same. This fact is not trivial, first because 
DD-ME2 has not been adjusted to nuclear matter data, 
but only to the experimental skin thickness of several fi- 
nite nuclei [l(| and second, the full isospin dependence 
is determined by the p-meson, whereas in DD-ME(5 it is 
distributed over the 6 and the p-meson. 

The reason for this good agreement can be understood 
from the upper panel of Fig. [5] where the different contri- 
butions to the symmetry energy are displayed, the kinetic 
part as well as those provided by the p and the <5-meson. 
We can see that the contributions of these mesons have 
opposite sign and thus a noticeable cancellation appears 
between them over the entire range of densities under 
consideration. Thus, it is conceivable (see Eq. (|40|) ) that 
if the (5-meson is not considered in the functional (as it 
is the case of DD-ME2) its contribution to the symme- 
try energy can be accounted for by the p-meson (with a 
reduced strength of the coupling constant, see Ref. [l(| 
and Table |TJ). 

The lower panel shows similar decompositions of the 
symmetry energy for other density functionals, such as 
NL3 [3 and DD-ME2 [Io|. The parameter set NL3 
(black) has no density dependence in the isovector chan- 
nel. Therefore the contribution of the p-meson is very 
stiff and proportional to the density. The parameter set 
DD-ME2 includes only one isovector meson, the p-meson 
and its contribution to the symmetry energy is very close 
the the sum of both the p- and the 5-meson for the set 
DD-ME<5 which compensate each other to a large extend. 
Small differences in these curves at densities above satu- 
ration density can be traced back to the different ansatz 
for the density dependence of the p-meson in these two 
parameter sets, the Typel-Wolter ansatz (|4T|) for DD- 
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FIG. 5: (Color online) Upper panel: The symmetry energy 
S2 (p) (full black) and its contributions as a function of the 
density, the kinetic contribution (dash dotted in blue), the 
contributions of the p-meson (dashed in red), the S- meson 
(dash dotted in green), and approximation (|37[) for the con- 
tribution of the 5-meson (dotted in yellow) . Lower panel: The 
symmetry energy S2(p) resulting from the parameter sets NL3 
(triangles), DD-ME2 (circles), and DD-ME5 (squares). The 
total value (solid line) is compared with the contribution of 
the p meson (dashed line) for NL3 and DD-ME2 and of the 
sum of p and 5 mesons for DD-ME5. 



B. Ground-state properties of finite nuclei 



As described in section IIII Al the experimental masses 
of 161 and the charge rms radii of 86 even-even spherical 
nuclei (see Table [Vj have been taken into account in the 
fitting procedure of the DD-ME5 functional. 

We display in Figs.|6]and[7]the difference between the- 
oretical results computed with the functionals DD-ME<5, 
DD-ME2 and BCP and experimental data. For DD-ME<5 
we obtain a rms deviation of 2.4 MeV for the binding 
energies and of 0.02 fm for the charge radii. These re- 
sults are close to the rms deviations 2.1 MeV and 0.02 
fm obtained with the DD-ME2 functional for the same 
set of data when pairing correlations are introduced by 
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the monopole force discussed at the end of Sect. IIII Bl 
It has to be emphasized, however, that using the den- 
sity functional DD-ME2 in connection with the pairing 
part of the finite range Gogny force D1S instead of the 
monopole force and taking into account spherical as well 
as deformed nuclei one has found rms deviations of 900 
kcV and 0.017 fm for the binding ener gies and charge 
radii of typical sets of 200 [l(| or 300 [l05| | even-even 
nuclei. 

The charge radii r c (defined in Eq. (|49I0 of Pb iso- 
topes and their isotope shifts have been a matter of de- 
tailed discussio n within the framework of mean field the- 
ories Il0(iiri0!> . In Fig. |S] we show the isotope shifts 
in a chain of Pb isotopes as a function of the neutron 
number N. The nucleus 208 Pb has been taken as the 



reference point: A 2 . (N) 



\(N) - ^(126). With a 



gradual addition of neutrons, the empirical charge radii 
of isotopes heavier than 208 Pb do not show the trend 
of the lighter isotopes and at the doubly mag ic nu- 
cleus 208 Pb one observes a pronounced kink llOj. Con- 
ventional non-relativistic Sky rme- and Gogny forces fail 
to reproduce this kink [106| . whereas all the relativis- 
ts model s are succes sful i n de scribing this kink prop- 
erly [To3. In Refs. [lOl [ToU this difference between 
the non-relativistic Skyrme functional and the rclativistic 
models has been traced back to the isospin dependence of 
the spin-orbit force. In conventional relativistic models 
it is determined by the p-meson vertex and it is relatively 
weak. In Fig. [S] we see that the parameter set DD-ME<5 
reproduces the kink in the isotope shifts rather well as all 
the other relativistic models do. The non-relativistic set 
BCP of Ref. [IH which has the same spin-orbit force as 
conventional Skyrme and Gogny functionals fails in this 
context. 



0.12 



■q 0.06 



I 

V-0.06 



-0.12 



iii! 

\ 
8 

"if" 
«□ * 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

BCP 

° DD-ME2 " 

* DD-ME5 _ 


~ x" n >P 

IP □ ~ 

© 

X © 

00 ©a 

© 

* .* 

„ 



1 1 1 1 


v 

an 

* 
□ 



50 



100 



150 



200 



FIG. 7: Color online) Difference between theoretical and ex- 
perimental charge radii as a function of the mass number. 
DD-ME5 results as compared with those of DD-ME2 [TT| and 
of BCP [T2| ■ The orange region corresponds to twice the fixed 
weight used in the fit (see Table IIII|) . 



Finally we show in Fig. [5] values for the neutron skin 
thickness Ar np = (r 2 )l/ 2 — {r 2 )^ 2 of a large set of nu- 
clei as a function of the relative neutron excess I = 
(N — Z)/ (N + Z) and compare the results obtained with 
the parameter set DD-ME<5 with those of the set DD- 
ME2 and with experimental values [9!| . Both theoretical 
calculations are in rather good agreement and within the 
range of the experimental error bars. 
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sponds to twice the fixed weight used in the fit (see Table 
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FIG. 8: (Color online) Isotope shifts for the chain of Pb iso- 
topes with respect to 208 Pb. Calculations with the relativistic 
model DD-ME5 and the non-relativistic functional BCP are 
compared. 
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C. Impact of the 5-meson on the spin-orbit 
splitting 

In this work, we have included the <5-meson in our the- 
oretical treatment of the nucleus motivated by micro- 
scopic calculations [H, 0, [65[ and by the importance of 
a scalar-isovector meson of the nucleon-nucleon poten- 
tials for descri bing the nucleon-nucleon scattering data 
in the vacuum [72| . In our investigation of the proper- 
ties of nuclear matter we have seen in Fig. [5] that the 
influence of the (5-mcson on the symmetry energy can be 
largely compensated by renormalizing the p-meson cou- 
pling constant in the DD-ME2 model. The same seems 
to be true also for the masses (Fig. |6j), radii (Fig. [7]) and 
skin thicknesses (Fig. |9]) in finite nuclei. Obviously this 
also applies for all the other successful covariant density 
functionals without the 5-meson degree of freedom. 

In order to g et a better understanding of these results 
we follow Ref. |lll| and eliminate the small components 
of the spinor ipi in the Dirac equation (J5J) . For the large 
components /i(r) we are left with a Schroedinger like 
equation 



1 



ap- 



<rp + V+\ f , 



2m + Ei + V- 
It contains the potentials 

v± = s s ± s° 



I ft ■ 



(54) 



(55) 



The potential V + « —50 MeV corresponds to the con- 
ventional potential in the corresponding non-rclativistic 
Schrodingcr equation. In theories containing the p and 
the (5-mesons it can be decomposed into an isoscalar and 
an isovector part 



Vl s (r) = g a v(r)+9 u u (r) 
Vl V {r)=gMr)+9 P Pl{r), 



(57) 
(58) 



where <r(r), uj (r), 5 3 (r), and p 3 (r) are the corresponding 
meson fields. In theories without (5-mcson the p-mcson 
vertex has to be renormalized and we find for the isovec- 
tor part a pure p-field 



(59) 



with a renormalized coupling g p . Since the field 5 3 (r) 
and Ps(r) have opposite sign the renormalized coupling 
g p has to be considerably smaller than the original g p , as 
it is seen in Tables Q] and [TTJ 
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FIG. 10: (Color online) The isoscalar part (left panel) and 
the isovector part (right panel) of the spin-orbit potential for 
the Sn-isotopes 100 Sn, 132 Sn, and 176 Sn. 
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The situation is different for the potential V- ~ —700 
MeV, which leads to a very large spin-orbit term. In 
spherical nuclei it has the strength U^°(r) for neutrons 
and protons (r = n,p) of the form 



1 



1 



IdVl 



2m 2m + VI r dr 



(60) 



Because of the non-linear connection between U so and V- 
thc situation is more complicated here and we obtain the 
isoscalar and isovector parts as 



U IS = - (U n + U p ) 

^ so 2 ^ so so J 



IV 



{Ko - u p so ) , 



(61) 
(62) 



Nonetheless these terms are dominated by V~(r) and the 
decomposition of this function with respect to isospin is 



V^(r)=g a a(r)-9^(r) 
Vl V (r)=g s 5 3 (r)-g pP °(r), 



(63) 
(64) 



Because of the opposite sign of the fields a and uj the 
corresponding isoscalar part of the spin-orbit potential is 
considerably enhanced with respect to the isoscalar part 
of the normal potential. This well known fact is also true 
for the isovector part. It is also considerably enhanced 
with respect to the isovector part of the normal potential. 
As a consequence there is an essential difference between 
a theory with and without a 5-meson, i.e. we expect an 
enhancement of the isospin-dependence of the spin-orbit 
potential in a theory with a i5-meson. 

To clarify these statements we show in Fig. [TTj] the 
isoscalar and the isovector part of the spin-orbit potential 
defined in Eqs. (pTTj) and (|6"2"|) for the three Sn-isotopcs 
100 Sn, 132 Sn. and 176 Sn. In order to have a clear evi- 
dence for the isospin dependence we neglect in this case 
the Coulomb potential. 

In the first panel ( 100 Sn) we have N = Z and because 
of neglecting the Coulomb potential the neutron and pro- 
ton densities are identical. As a consequence the isovec- 
tor meson fields 5(r) and p(r), as well as the isovector 
part of the spin-orbit potential vanish identically. This 
is true for both models DD-MEJ and DD-ME2. Because 
of the different coupling constants in these models the 
isoscalar parts of the spin-orbit potentials are slightly 
different, but large in both cases and peaked at the sur- 
face. In the second panel ( 132 Sn) we have a considerable 
neutron excess. The isoscalar part did not change very 
much. Apart from the fact that the larger mass number 
A = 132 produces a shift of the surface and the maximum 
of the spin-orbit potential to larger r-values, both mod- 
els show similar results. The situation is very different 
for the isovector part. For the DD-ME2 model without 
(5-meson the p-field is not vanishing but relatively small 
as compared to the isoscalar fields (please consider the 
change of the scale in the right hand side of this figure). 
The effective coupling for the p-exchange value of g 2 p is 



considerably smaller than the other two couplings for the 
other two mesons <r or w. In addition the source of the p 
field is the difference between the neutron and the proton 
density. This difference is even for N=82 not extremely 
large. On the other side for the DD-ME5 model with a 
6 and a p meson the isovector part of the spin-orbit field 
is considerably enhanced with respect to the isovector 
part of the DD-ME2 model. However, since this is rela- 
tively small, in total the isovector part of the spin-orbit 
potential is still an order of magnitude smaller than the 
isoscalar part. Therefore even for effects, which depend 
on the spin-orbit potential we do not expect essential dif- 
ferences between a model with and without a (5-meson. 
This is even true for cases with extreme neutron excess 
as in the nucleus 176 Sn in the lowest panel of Fig. [TU1 

In Fig. [TT]wc show the spin-orbit splitting for neutron 
orbitals in a chain of Sn-isotopes, starting at the N = Z 
nucleus 100 Sn. In contrast to Fig. [10] the Coulomb inter- 
action is included here. Results for the two parameter 
sets DD-ME2 (without 5-meson) and DD-ME5 (with 5- 
meson) are compared. First we find that the difference 
of these two models is rather small for the 2p- and the 
2d-orbits. Because of the low £-values these splittings 
are relatively small and the corresponding wave functions 
are not so surface peaked. For the If - and the lg-orbits 
with large £-values the splitting is relatively large and 
we find a considerable difference between the DD-ME2 
and the DD-MEc> model. This difference is, however, 
connected with the isoscalar part of the spin-orbit po- 
tential, because it occurs already in the nucleus 100 Sn 
which has, apart from a small violation of isospin due to 
the Coulomb force, practically no isovector part. 

With increasing neutron number the spin-orbit split- 
ting in these high I orbitals decreases. This ha s al- 
read y bee n observed in earlier investigations in Ref. jl!2| 
and [113j, where it has been explained by the increasing 
neutron skin and the increasing neutron diffuseness lead- 
ing to a reduced derivative in the spin-orbit potential. 

Finally we observe in Fig. QT] that the difference in the 
neutron spin-orbit splittings calculated with DD-ME2 
and DD-ME<5 decreases slightly with increasing neutron 
number. This effect has its origin in the increasing isovec- 
tor part of the neutron spin-orbit splitting for the param- 
eter set DD-ME<5: [7™ = (C/ S 7 s + C// y )/2. However, it is 
relatively small because the isovector part itself is small 
as compared to the isoscalar part (see the scales on the 
right panels in Fig. [TTj] Of course, for the protons with 
^so = (Uio — Ulo)/2 this difference increases with in- 
creasing neutron number (not shown in Fig. 1111) . 



V. DISCUSSION AND CONCLUSIONS 

On the way to a more microscopic derivation of rel- 
ativistic nuclear energy density functionals, we started 
with Brucckner calculations [l9( for symmetric nuclear 
matter and pure neutron matter and with Dirac Brueck- 
ner calculations fo5j for pure neutron matter. We tried to 
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FIG. 11: (Color online) Spin-orbit splitting in Sn-isotopes 
for various nl levels. Results for the parameter set DD-ME5 
(green crosses) are compared with those of the parameter set 
DD-ME2 (black circles) 



use this microscopic information as far as possible for the 
adjustment of a new covariant density functional based 
on density dependent meson exchange. Since it is well 
known that, at present, all attempts to derive the func- 
tionals directly from bare forces do not reach the required 
accuracy for nuclear structure applications we added ex- 
perimental data in finite nuclei, such as binding; energies 
and charge radii for the fit. In contrast to Ref. [16| , where 
a similar idea has been applied to a relativistic point cou- 
pling model, we took into account in this work the fact 
that in Dirac Brueckncr calculations the resulting scalar 
self energies show a strong isovector part by including a 
(5-meson which is usually neglected in relativistic meson 
exchange models. 

This investigation is in some sense an extension of ear- 
lier non-relativistic work in Refs. [Hj]. We want to point 
out, however, that in the present work we were forced to 
leave, in part, the strategy o f Re f. (l2j . In that paper 
and in earlier work of Fayans [l!4| the strict Kohn-Sham 
strategy was followed as in Coulombic systems. Namely 
the bulk part was exclusively determined from previous 
microscopic Brueckner calculations [lj| and thus fixed 
once and for all. Then a phenomenological finite range 
contribution and a spin-orbit term was added to the func- 
tional to account for properties of finite nuclei in ad- 
justing four parameters. Since in the relativistic case 
the spin-orbit is fixed already from the nuclear matter 
calculations (37|, if all the parameters of the relativis- 
tic approach which survive in the infinite matter limit 
were adjusted to microscopic infinite matter results, one 
would essentially remain with only one adjustable pa- 
rameter, i.e. the tr-mass, which can serve for the ad- 
justment to properties of finite nuclei. The u- and the 
p-mass cannot be fixed independently from present data 
and therefore they are kept at their experimental values 
in the vacuum. It turns, however, out that only one pa- 
rameter is not enough to reach the required accuracy for 



nuclear masses and radii. We, therefore, had to adjust 
parameters simultaneously to microscopic bulk proper- 
ties and experimental finite nuclei data, thus departing 
from the strict Kohn-Sham strategy. Nevertheless, we 
ensured that in our model the density dependence of the 
meson exchange couplings is completely determined by 
the microscopic calculations of infinite matter and only 
four remaining parameters are adjusted to experimental 
data in finite nuclei. 

On the other side we have to keep in mind that the nu- 
clear many-body problem is much more complicated than 
Coulombic systems. It is based on QCD, a relativistic 
theory, where spin-degrees of freedom play an essential 
role, not only in the spin-orbit term. Relativistic mod- 
els provide a consistent treatment of the spin degrees of 
freedom and velocity dependent terms, they include the 
complicated interplay between the large Lorentz scalar 
and vector self-energies induced on the QCD level by the 
in mediu m ch anges of the scalar and vector quark con- 
densates [l!5j . In particular they include nuclear mag- 
netism, i.e. the nuclear currents induced by the spatial 
parts of the vector self energies or the time-odd com- 
ponents of the nuclear density density functional. The 
Kohn-Sham strategy applied in Ref. [12| has, at present, 
no possibility to derive this part of the functional. 

As a result of our investigations we have derived a func- 
tional DD-ME£ with similar properties as the very suc- 
cessful functional DD-ME2. In contrast to that model, 
DD-ME(5 is based to a large extent on microscopic cal- 
culations. Only four parameters had to be adjusted to 
finite nuclei. It turns out that the inclusion of the 5- 
meson does not improve the accuracy of the properties 
of finite nuclei such as masses and radii. Therefore the 
corresponding vertex and its density dependence is com- 
pletely determined by nuclear matter data such as the 
isovector part of the effective Dirac mass. It is neverthe- 
less much more physical and notably the mass splitting 
of neutrons and protons is now correctly incorporated. It 
has, moreover, an influence on the behavior of the equa- 
tion of state at higher densities and we find in this region 
a much better agrement with experimental data derived 
from heavy ion reactions [H| than the earlier parameter 
set DD-ME2 not including the (5-meson. We therefore 
can hope that the new parameter set derived in this in- 
vestigation is more reliable for applications of relativistic 
density functional theory to neutron stars. 
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TABLE V: DD-ME5 results for the binding energies and 
charge radii used in the fit 



